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Preface 


This  report  presents  the  results  of  a  numerical  model  investigation  used  to 
predict  average  salinity  changes  that  will  occur  in  the  Lake  Pontchartrain  Basin  as 
a  result  of  varying  levels  of  closure  of  the  Mississippi  River-Gulf  Outlet  below 
Lake  Borgne,  LA. 

This  investigation  was  conducted  from  January  2000  through  August  2000  at 
the  U.S.  Army  Engineer  Research  and  Development  Center  (ERDC),  Vicksburg, 
MS,  by  Dr.  R.  C.  Berger,  Mr.  A.  R.  Carrillo,  and  Ms.  M.  S.  Samiff  of  the  Coastal 
andHydraulics  Laboratory  (CHL),  and  Mr.  B.  J.  Thibodeaux  of  the  U.S.  Army 
Engineer  District,  New  Orleans  (LMN).  Funding  was  provided  by  the  New 
Orleans  District. 

The  work  was  performed  under  the  general  direction  of  Mr.  Thotnas  W. 
Richardson,  Acting  Director,  CHL,  and  Dr.  R.  T.  McAdory,  Chief,  Tidal 
Hydraulics  Branch. 

At  the  time  of  publication  of  this  report,  Dr.  James  R.  Houston  was  Director 
of  ERDC,  and  COL  John  W.  Morris  IB,  EN,  was  Commander. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication, 
or  promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 
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1  Introduction 


Background 

The  Mississippi  River-Gulf  Outlet  (MRGO)  consists  of  a  ship  channel 
36  feet  deep  and  500  feet  wide,  extending  approximately  76  miles  from  the 
junction  of  the  Inner  Harbor  Navigation  Canal  and  the  Gulf  Intracoastal 
Waterway  in  New  Orleans,  LA,  to  the  -38  mlw  (mean  low  water)  foot  contour  in 
the  Gulf  of  Mexico.  The  purpose  of  the  MRGO  is  to  provide  a  deep  draft 
channel  to  the  Port  of  New  Orleans  Inner  Harbor  facilities. 

Since  the  MRGO's  completion  in  January  1968,  saltwater  flux  from  the 
MRGO  through  direct  connections  to  Lake  Borgne  and  the  Gulf  Intracoastal 
Waterway  has  contributed  to  an  increase  in  the  salinity  concentration  of  the  lakes 
and  Biloxi  Marshes.  Attempts  by  the  U.S.  Army  Engineer  District,  New  Orleans 
(LMN)  to  mitigate  this  salinity  increase  through  the  operation  and  maintenance 
program  have  met  with  little  or  no  success. 

Therefore,  the  New  Orleans  District  is  participating  in  an  Environmental 
Protection  Agency  sponsored  study,  “Mississippi  River  -  Gulf  Outlet,  LA, 
Re-evaluation  Study,”  focusing  on  the  deep  draft  navigation,  environmental,  and 
flood  control  aspects  of  the  project  and  the  need  for  continued  maintenance.  The 
purpose  of  the  re-evaluation  study  is  to  determine  if  the  existing  ch^el  should 
be  modified  and  to  provide  the  advisability  of  continuing  its  operation.  In 
support  of  the  effort,  the  U.S.  Army  Engineer  Research  and  Development  Center, 
Waterways  Experiment  Station  (ERDC-WES)  Coastal  Hydraulics  Laboratory 
(CHL)  was  tasked  to  model  several  plans  for  partial  blockage  of  the  MRGO 
south  of  Lake  Borgne,  LA. 


Objective 

The  objective  of  tiie  work  presented  herein  is  to  predict  the  average  salinity 
changes  that  will  occur  in  the  Mississippi  and  Louisiana  estuaries  of  the  Lake 
Pontchartrain  Basin  (Figure  1),  particularly  that  part  known  as  the  Biloxi 
Marshes,  resulting  from  partial  closure  of  the  MRGO  below  Lake  Borgne. 

The  purpose  of  this  report  is  to  present  the  results  of  a  numerical  model 
investigation  addressing  this  objective.  The  modeled  region  includes  Lake 
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Pontchartrain,  Lake  Borgne,  the  MRGO,  the  hiner  Harbor  Navigation  Channel 
(IHNC),  the  Gulf  Intracoastal  Waterway  (GIWW),  the  Rigolets,  the  Chef 
Menteur,  segments  of  Mississippi  Sound,  Chandeleur  Sound,  Breton  Sound,  and 
at  least  a  portion  of  the  Biloxi  Marshes,  as  shown  in  Figure  1 . 

Based  on  the  results  of  the  initial  study,  additional  investigation  was 
requested  to  evaluate  the  length  of  time  necessary  for  the  system’s  salinity  to 
return  to  normal  observed  levels  after  an  unusual  weather  event.  In  this  case,  the 
surge  from  Hurricane  Juan  is  used  to  constitute  this  event.  The  results  from  this 
additional  investigation  can  be  found  in  Appendix  B. 
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Lake  Pontchartrain  Basin 


The  Lake  Pontchartrain  Basin  consists  of  Lake  Maurepas,  Lake 
Pontchartrain,  Lake  Borgne,  the  Biloxi  Marshes,  and  Chandeleur  Sound,  plus 
associated  marshlands  and  waterways.  The  Basin  is  described  in  detail  by  the 
New  Orleans  District  (USAED,  New  Orleans,  1984, 1990),  Pankow  et  al.  (1989), 
and  the  U.S.  v^rmy  Corps  of  Engineers  Committee  on  Tidal  Hydraulics  (CTH) 
(1995).  A  summary  of  pertinent  factors  will  be  provided  here. 


Hydrology 

The  largest  tributary  to  the  area  is  the  Pearl  River  with  a  mean  annual  flow  of 
about  10,000  cfs.  The  Pearl  River  discharges  into  Lake  Borgne  near  the  mouth 
of  the  Rigolets,  one  of  three  tidal  waterways  out  of  Lake  Pontchartrain.  Several 
smaller  rivers,  the  largest  of  which  are  the  Amite,  Tickfaw,  Tangipahoa,  and 
Tchefuncta  Rivers,  flow  into  Lake  Pontchartrain  and  Lake  Maurepas.  The  annual 
freshwater  flow  into  Lake  Pontchartrain  averages  about  3,800  cfs. 


Hydrodynamics 

Tides  in  the  basin  are  principally  diurnal,  with  mean  ranges  of  0.3  ft 
(Lake  Maurepas)  to  0.5  ft  (Lake  Pontchartrain)  to  1.4  ft  (Chandeleur  Sound). 
Sustained  winds  can  raise  or  lower  peak  astronomical  tide  levels  by  several  feet 
for  short  periods.  Mean  water  levels  are  affected  by  winds,  freshwater  runoff, 
and  seasonal  trends  in  the  Gulf  of  Mexico. 

Figure  2  illustrates  typical  seasonal  variability  in  mean  water  level.  It 
displays  monthly  mean  water  levels  for  1940-1950  at  Eugene  Island  at  the  Gulf 
entrance  to  Atchafalaya  Bay  (USAED,  New  Orleans,  1982)  which  are 
representative  of  Gulf  levels  throughout  the  area  in  that  mean  water  levels  are 
lower  in  winter  and  midsummer  and  higher  in  spring  and  fall. 

Currents  and  circulation  are  controlled  by  tides,  winds,  freshwater 
discharges,  and  gulf  currents.  Flows  in  the  ^^GO  are  also  affected  to  some 
extent  by  density  currents  (Donnell  and  Letter  1991). 
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Figure  2.  Intra-annual  mean  water  level  variation  at  Eugene  Island,  1940-50 
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Lake  Pontchartrain  Historical  Saiinity 

Historical  salinity  data  for  five  stations  in  the  Pontchartrain  Basin  were 
analyzed.  These  stations  are  Pass  Manchac  near  Pontchatoula,  Lake  Pontchar¬ 
train  at  Little  Woods,  Chef  Pass  near  Lake  Borgne,  Lake  Pontchartrain  at  North 
Shore  and  Bayou  LaLoutre  at  Alluvial  City.  The  data  are  evaluated  as  pre-  and 
post-MRGO  periods  to  identify  changes  to  the  salinity  regime  due  to  the  MRGO. 
The  bulk  of  the  project  was  completed  in  July  1963,  but  a  plug  at  Paris  Road  was 
not  removed  until  1968.  Nevertheless,  the  condition  between  1963  and  1968  did 
allow  intrusion  up  the  MRGO  to  Lake  Borgne. 


The  data  indicate  that  the  lowest  salinities  are  generally  in  the  late  spring  and 
highest  in  the  summer  and  late  fall.  This  reflects  seasonal  variations  in 
freshwater  inflows  from  the  major  rivers  and  streams  in  the  basin.  The  salinities 
of  Lake  Borgne  and  Lake  Pontchartrain  normally  range  from  fresh  to  brackish. 
During  periods  of  extreme  low  flows.  Lake  Maurepas  can  become  brackish. 

Salinities  in  Lake  Borgne  generally  range  from  2  to  1 5  ppt  and  are  strongly 
influenced  by  Pearl  River  discharges  and  inflow  from  the  Rigolets  and  Chef  Pass 
Higher  salinity  water  from  the  MRGO  enters  Lake  Borgne  through  breaks  in  the 
marshes  between  the  two  water  bodies. 
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Analyses  of  salinity  data  indicate  that  the  most  notable  increase  in  monthly 
average  salinity  occurred  after  1963.  Mean  monthly  salinity  increased  for  all 
months  for  the  period  subsequent  to  1963.  This  increase  can  be  attributed  partly 
to  the  partial  completion  of  the  MRGO  in  1963  which  provided  a  major  access 
for  saline  water  to  enter  Lake  Maurepas,  Lake  Pontchartrain,  and  Lake  Borgne. 

Monthly  summaries  of  salinity  for  pre-  and  post-MRGO  indicate  that  salinity 
has  increased  on  the  average  by  the  following  amounts: 

a.  1 .3  ppt  at  Lake  Pontchartrain,  North  Shore. 

b.  1.8  ppt  at  Lake  Pontchartrain,  Little  Woods. 

c.  0.4  ppt  at  Pass  Manchac  near  Pontchatoula. 

d.  2.4  ppt  at  Chef  Pass  near  Lake  Borgne. 

e.  4.3  ppt  at  Bayou  LaLoutre,  Alluvial  City. 

The  salinity  regime  in  ftie  area  has  stabilized  and  no  sigmficant  increase  in 
average  annual  salinity  is  projected  in  the  foreseeable  futxire  for  Lake  Maurepas 
and  Lake  Pontchartrain.  Salinity  is  expected  to  increase  in  the  Lake  Borgne  area 
and  surrounding  marshes  due  to  land  loss  in  the  area. 

Table  1  shows  mean  monthly  pre-  and  post-MRGO  salinity  for  the  period 
1951  to  1963  and  1963  to  1977. 


Table  1 

Pre-  and  Post-MRGO  Salinity  ppt 
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Numerical  Model 


TABS-MDS,  the  numerical  model  used  for  this  investigation,  is  a  three- 
dimensional  finite  element  code  originally  developed  by  Dr.  Ian  King  of 
Resource  Management  Associates,  and  modified  at  ERDC.  It  models  three- 
dimensional  hydrodynamics  and  salt  transport  accounting  for  unsteady  river 
inflows,  tides,  wind  effects,  and  density-driven  circulation.  It  has  been  widely 
used  by  ERDC  to  model  three-dimensional  hydrodynamics  and  salinity  at 
numerous  locations,  including  Galveston  Bay,  TX  (Berger  et  al.  1995).  The 
model  is  described  more  fully  in  Appendix  A. 


Computational  Mesh 


Figure  3  illustrates  the  planform  view  of  the  computational  mesh  used.  The 
mesh  is  three-dimensional  everywhere  except  near  the  gulfward  boundary  and  in 
Lake  Maurepas.  The  3-D  grid  contains  16,932  nodes  and  6,701  elements 
distributed  horizontally  and  vertically.  The  grid  is  refined  along  the  MRGO,  and 
its  connections  to  the  Gulf  of  Mexico,  Lake  Borgne,  and  Lake  Pontchartrain.  The 
computational  mesh  provides  direct  connections  between  MRGO  and  Lake 
Borgne  at  Shell  Beach  and  Martello  Castle.  Another  connection  exists  between 
Lake  Borgne  and  the  GIWW  near  Bayou  Gentilly  (southeast  of  Chef  Menteur). 
These  connections  are  sized  to  approximate  not  only  the  connections  at  those 
locations,  but  also  nearby  smaller  connections,  and  thus  represent  an  aggregation 
of  several  smaller  waterways. 


Model  Validation 

Validation  for  the  Lake  Pontchartrain  Basin  model  was  performed  in  a  prior 
investigation.  Details  of  the  validation  can  be  found  in  McAnally  and  Berger 
(1997).  ^  ^ 
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Figure  3.  Computational  mesh  and  region  of  closure 


Experimental  Conditions 

Typical  modeling  conditions  for  the  base  and  plan  experiments  are  based  on 
statistical  measures.  Since  most  of  this  work  is  based  on  prior  investigations, 
only  a  brief  discussion  of  the  specifics  of  the  experimental  conditions  are 
addressed  in  the  following  sections.  For  a  more  detailed  discussion  see 
McAnally  and  Berger  (1997). 


Boundary  Conditions 

Riveiflows  for  base  and  plan  experiments 

Base  and  plan  experiments  were  conducted  for  typical  conditions,  using  the 
50  percent  exceedance  flows  given  in  Table  2.  These  have  been  corrected  for 
ungauged  areas.  A  complete  description  of  the  development  of  the  flows  used 
for  the  model  (Table  2)  is  given  in  McAnally  and  Berger  (1997). 


Tides 

Tides  at  the  model's  gulfward  boundary  were  synthesized  for  the  year  1982 
from  the  tidal  constituents  as  given  by  Outlaw  (1982)  and  shown  in  Table  3.  The 
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Table  2 
Stream  1 
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Table  3 
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S2 
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^5.0 

N2 
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0.02 

209.4 

0.02 

223.0 

0.01 

183.8 

mean  water  level  set  to  zero  referred  to  the  National  Geodetic  Vertical  Datum 
(NGVD). 


Wind 

The  wind  data  used  were  obtained  from  the  U.S.  Air  Force  Environmental 
Technical  Applications  Center  in  Asheville,  NC.  These  data  were  the  hourly 
surface  winds  at  the  New  Orleans  International  Airport  for  the  calendar  year 
1 982,  and  were  used  for  all  base  and  plan  experiments. 


Initial  Conditions 

In  these  experiments  an  initial  salinity  field,  currents,  and  water  elevations, 
were  obtained  from  previous  simulations  of  the  region  (McAnally  and  Berger  ’ 
1997),  ^d  used  as  initial  conditions  for  all  base  and  plan  experiments.  Since  the 
simulations  were  started  in  January,  the  initial  conditions  had  sufficient 
relaxation  time  before  the  first  period  of  interest  in  April. 
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Base  and  Plans 

The  base  condition  (maximum  deptii  in  the  MRGO  cross  section  is  47  ft)  and 
four  levels  of  closure  were  modeled.  Closure  consisted  of  a  sill  across  the 
LaLoutre  Ridge  in  the  MRGO,  and  the  region  is  labeled  in  Figure  3 .  All 
employed  the  boundary  and  initial  conditions  described  previously.  The 
experiments  were  as  follows; 

a.  Base:  Maximum  depth  in  the  cross  section  is  47  ft. 

b.  20-ft  Plan:  Closed  to  within  a  20-ft  depth. 

c.  1 6-ft  Plan:  Closed  to  within  a  16-ft  depth. 

d.  12-ft  Plan;  Closed  to  within  a  12-ft  depth. 

e.  Complete  closure. 

The  numerical  model  calculated  water-surface  elevations,  current  velocities 
(three-dimensional  components),  and  salinities  at  each  node  every  60  minutes  for 
the  10-month  period  of  simulation  between  January  1  and  October  3 1 .  Those 
data  were  processed  to  provide  average  monthly  salinity  contour  plots  for  base 
and  each  plan  for  April,  May,  September,  and  October  (Plates  1-4). 


Chapters  Approach 


4  Results  and  Discussion 


The  interest  of  this  study  is  the  effect  of  the  closure  of  MRGO  on  salinities  in 
Lake  Pontchartrain,  Lake  Borgne,  and  Biloxi  Marsh.  The  results  are  contained  in 
Plates  1-4.  These  contain  the  month  average  isohalines  for  the  bottom  depth  for 
each  plan.  The  isohalines  shown  for  each  plan  represent  the  change  from  the 
base  conditions.  The  base  conditions  are  the  month  average  salinity.  The  plan 
isohalines  are  then  changes,  where  a  negative  sign  (-)  indicates  that  the  closure 
reduced  the  salinity  and  a  positive  sign  (+)  indicates  an  increase.  The  month 
averages  given  are  for  April,  May,  September,  and  October.  Tables  4-7  show  the 
values  for  specific  station  locations  (approximate).  The  spring  months  are 
representative  of  the  low  salinity  period  of  the  year,  and  the  autumn  months 
indicate  the  maximum  salinity  period. 

The  base  isohalines  are  found  on  part  (a)  of  Plates  1, 2, 3,  and  4.  The  spring 
months  show  salinity  values  in  Lake  Pontchartrain  that  are  between  4  and  6  ppt. 
In  Lake  Borgne  the  spring  month  average  salinity  is  generally  shown  to  be 
between  8  and  16  ppt.  The  fall  month-average  salinity  values  in  Lake 
Pontchartrain  are  generally  between  4  and  8  ppt.  In  Lake  Borgne  in  the  same 
period,  the  salinity  values  are  shown  to  be  roughly  between  12  and  20  ppt. 

An  obvious  remark  that  can  be  made  about  the  salinity  for  each  plan  is  that 
the  largest  changes  result  from  the  complete  closure.  Even  the  12-ft  plan,  which 
involves  leaving  a  MRGO  blockage  for  depths  below  12  ft  mean  low  water 
(mlw),  is  much  less  dramatic  than  the  complete  closure.  The  major  result  with 
any  closure  is  that  Lake  Borgne  and  Lake  Pontchartrain  salinities  decrease. 

There  is  veiy  little  change  in  salinity  along  the  eastern  edge  of  Lake  Borgne.  The 
largest  salinity  decreases  occur  along  the  western  portion  of  Lake  Borgne  and  the 
southeastern  shore  of  Lake  Pontchartrain.  These  are  a  result  of  the  two  major 
connections  to  MRGO  from  Lake  Borgne  and  the  Inner  Harbor  Navigation  Canal 
that  connect  Lake  Pontchartrain  to  MRGO. 


The  average  salinity  reduction  for  any  of  the  partial  closure  plans  was  0.3  ppt 
or  less  in  Lake  Pontchartrain  and  Lake  Borgne.  For  the  complete  closme  the 
spring  months  showed  general  salinity  decreases  compared  to  base  of  between 
0  and  3.4  ppt  in  Lake  Borgne  and  0.2  to  2.4  ppt  in  Lake  Pontchartrain.  The 
autumn  months  showed  general  decrease  of  0.2  to  3.4  ppt  in  Lake  Borgne  and 
between  1.0  to  3.4  in  Lake  Pontchartrain.  The  Biloxi  Marsh  area  of  eastern 
Lake  Borgne  showed  very  little  change. 
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Table  4 

Modeled  April  Monthly  Average  Salinity,  Base 
Chanoes  from  Base,  ppt 

Condition 

s  and  Pian 

Location 

Base 

20-ft  Plan 

16-ft  Plan 

12-ft  Plan 

Closure  Plan 

Alluvial  City 

16.5 

-0.5 

-0.5 

-0.6 

-6.0 

Chef  Pass 

8.4 

-0.2 

-0.2 

-0.6 

-1.7 

Fenier 

4.6 

-0.0 

-0.1 

-0.1 

-0.6 

Little  Woods 

6.9 

-0.2 

-0.2 

-0.3 

-1.6 

Martello  Castle 

15.1 

-0.7 

-0.7 

-0.7 

-6.6 

North  Shore 

5.4 

-0.1 

-0.1 

-0.2 

-0.9 

Pass  Manchac 

0.7 

-0.0 

-0.0 

-0.0 

-0.1 

Pointe  Aux  Marchettes 

13.9  1 

-0.1 

-0.0 

0.0 

-0.5 

Table  5  ^ 

Modeled  May  Monthly  Average  Salinity,  Base  Conditiont 

Chanoes  from  Base,  ppt 

and  Plan 

Location 

Base 

20-ft  Plan 

16-ft  Plan 

12-ft  Plan 

Closure  Plan 

Alluvial  City 

l6?i 

-0.4 

-0.5 

-5.7 

Chef  Pass 

8.9 

-0.2 

-0.2 

-0.3 

-2.2 

Fenier 

4.7 

-0.1 

-0.1 

-0.1 

-0.8 

Little  Woods 

6.2 

-0.3 

-0.3 

-0.4 

-2.1 

Martello  Castle 

15.1 

-0.6 

-0.7 

-0.8 

-6.6 

North  Shore 

5.7 

-0.2 

-0.2 

-0.2 

-1.2 

Pass  Manchac 

0.6 

-0.0  ^ 

-0.0 

-0.0  ^ 

-0.1 

Pointe  Aux  Marchettes 

14.3 

-0.1 

-0.1 

-0.0 

-0.8 _ 

Table  6 

Modeled  September  Monthly  Average  Salinity,  Base  Conditions 
and  Plan  Chanoes  from  Base,  ppt 

Location 

Base 

20-ft  Plan  1  16-ft  Plan  |  12-ft  Plan  Closure  Plan  J 

Alluvial  City 

17.9 

-0.4 

-0.4 

-0.6 

-5.1 

Chef  Pass 

10.5 

-0.2 

-0.2 

-0.2 

-1.6 

Fenier 

4.9 

-0.1 

-0.2 

-0.2 

-1.3 

Little  Woods 

7.1 

-0.2 

-0.3 

-0.3 

-2.3 

Martello  Castle 

16.7 

-0.4 

-0.5 

-0.7 

-5.4 

North  Shore 

6.9 

-0.2 

-0.2 

-0.2 

-1.3 

Pass  Manchac 

1.0 

-0.0 

-0.0 

-0.3  ^ 

-0.3 

Pointe  Aux  Marchettes 

15.8 

-0.1 

-0.1 

-0.0 

•0.8 

Table  7 


Modeled  October  Monthly  Average  Salinity,  Base  Conditions  and 


Location 

Base 

20-ft  Plan  1  16-ft  Plan  1  12-ft  Plan  I  Closure  Plan  Jj 

Alluvial  City 

20.2 

-0.5 

-0.6 

-0.7 

-6.6 

Chef  Pass 

11.7 

-0.2 

-0.2 

-0.3 

-1.9 

Fenier 

5.4 

-0.2 

-0.2 

-0.2 

-1.6 

Little  Woods 

8.1 

-0.3 

-0.3 

-0.4 

-3.1 

Martello  Castle 

19.3 

-0.5 

-0.6 

-0.8 

-7.2 

North  Shore  1 

7.4 

-0.2 

-0.2 

-0.2  n 

-1.4 

Pass  Manchac 

1.1 

-0.0 

-0.0 

-0.1 

-0.4 

Pointe  Aux  Marchettes 

17.3 

-0.1 

-0.1 

-0.1 

-1.1 
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It  is  worth  noting  that  the  change  in  salinity  accompanying  the  complete 
closure  of  MRGO  is  similar  to  that  expected  from  the  historical  salinity  data. 
Table  8  shows  the  measured  and  modeled  changes  in  salinity  from  pre-  to  post- 
MRGO  for  five  stations  in  the  Pontchartrain  Basin.  The  measured  values  are 
based  on  the  mean  monthly  salinity  for  the  period  1951  to  1963  (pre-MRGO)  and 
the  period  1963  to  1977  (post-MRGO).  The  modeled  values  were  computed 
using  median  flows  rates,  astronomical  tides  with  limited  tidal  periods,  and 
1982  wind  data.  Given  these  differences,  the  model  results  are  quite  reasonable. 


Table  8 
Measuret 
MRGO,  p 

and  Mo 
Pt 

deled  Changes  in  Salinity  from  Pre-  to  Post- 

Month 

Pass 

Manchac 

North 

Shore 

Little 

Woods 

Chef 

Pass 

April 

Measured 

0.5 

0.7 

1.6 

1.8 

Modeled 

0.1 

0.9 

1.6 

1.7 

6.0 

May 

Measured 

0.1 

0.3 

1.6 

1.4 

0.7 

Modeled 

0.1 

0.1 

2.1 

2.2 

5.7 

September 

Measured 

0.3 

2.1 

1.7 

2.6 

4.7 

Modeled 

0.3 

1.3 

2.3 

1.6 

5.1 

October 

Measured 

0.4 

2.6 

1.9 

3.2 

6.2 

Modeled 

0.4 

1.4 

3.1 

1.9 

6.6 

Chapter  4  Results  and  Discussion 


5  Conclusions 


This  investigation  is  concerned  with  various  stages  of  closure  of  the  MRGO 
channel  from  the  Gulf  of  Mexico  to  the  city  of  New  Orleans.  IDstorical  records 
indicate  that  when  the  channel  was  built  there  was  an  increase  in  salinity  in  Lake 
Pontchartrain  and  Lake  Borgne.  This  model  study  used  a  sill  near  the  connection 
of  MRGO  to  the  Gulf  of  Mexico  with  an  elevation  of -20,  -16,  and  -12  ft  mlw  to 
provide  partial  closure  of  the  channel  and  investigate  the  restoration  of  the 
historical  salinity  regimen.  The  model  study  also  includes  the  base  conditions  of 
a  fully  open  channel  and  complete  closure  of  the  MRGO. 

The  monthly  average  salinities  from  the  model  indicate  that  all  of  the  closure 
plans  reduced  the  salinity  of  Lake  Borgne  and  Lake  Pontchartrain.  None  made 
large  changes  in  tibie  salinity  in  the  Biloxi  Marshes.  The  change  is  most 
pronounced  in  Lake  Pontchartrain  near  the  Inner  Harbor  Navigation  Canal, 
which  is  the  connection  to  MRGO.  The  eastern  portion  of  Lake  Pontchartrain 
had  a  more  significant  reduction  in  salinity  than  did  the  western  portion.  In  Lake 
Borgne,  salinities  in  the  areas  near  the  connection  to  MRGO  were  most  strongly 
reduced. 

The  partial  closures  had  a  very  small  effect.  Even  the  MRGO  closure  plan 
that  had  the  sill  elevation  of  only  -12ft  mlw  produced  reductions  that  were  quite 
small  compared  to  complete  closure. 

The  final  point  is  that  the  model  results  for  complete  closure  compared  to 
base  conditions  were  consistent  with  the  field  historical  salinity  changes  observed 
from  before  to  after  MRGO  was  built. 
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Appendix  A 

The  Hydrodynamic  Code 


The  complexity  of  this  estuary  requires  a  numerical  model  that  relies  upon 
an  unstructured  computational  mesh.  The  code  chosen  is  the  Galerkin-based 
finite  element  model  TABS-MDS,  which  is  a  U.S.  Army  Engineer  Research  and 
Development  Center  (ERDC)  adaptation  of  the  RMA-10  code  developed  by 
King  (1993).*  This  code  computes  time-varying  open-  channel  flow  and 
salinity/temperature  transport  in  1 , 2,  and  3  dimensions.  It  invokes  the 
hydrostatic  pressure  and  mild  slope  assumption.  Vertical  turbulence  is  supplied 
using  a  Mellor-Yamada  Level  H  (Mellor  and  Yamada  1982)  k-i  approach  modi¬ 
fied  for  stratification  by  the  method  of  Henderson-Sellers  (1984).  The  salin- 
ity/density  relationship  is  based  upon  Pritchard  (1982). 

The  fixll  three-dimensional  equations  are  reduced  to  a  set  of  two  momentum 
equations,  an  integrated  continuity  equation,  a  convection-diffusion  equation, 
and  an  equation  of  state.  The  simplification  is  a  result  of  the  hydrostatic 
pressure  approximation. 
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Dt  dx 


r  =0 


Oy 


-ry=o 


8h 


+ 


dt 


dC  da^  di; 

‘^~dx'^‘’dx  dy  dy  dy 


\dz=  0 


Dt 


dy  dy  )  dz 


ds\_^ 

D,—  =0 

uZ  y 


P=F(s) 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 


‘  All  references  cited  in  this  appendix  are  listed  in  the  References  section  at  the  end  of  the  main  text. 
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A1 


Elevation-related  terms  are  defined  in  Figure  A1 . 


Figure  A1.  Definitions  for  elevation  terms 
where 


1 

Exx  ^ 

Eyx  ^ 

dx 

dx 

{  Exy  ^ 

}■  ;  cr  y  = 

{  Eyy  ^  } 

dy 

Exz  ^ 

Eyz  ^ 

1  dz 

J 

1  dz  J 

and 

p  =  density 

u,v,w  =  x,y,z  velocity  components 
t  =  time 
P  =  pressure 

=  2a)sin(p) 

o)  =  rate  of  angular  rotation  of  the  earth 
<p  =  local  latitude 


A2 
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g  =  gravitational  acceleration 

(j  =  Chezy  or  Manning  friction  formulation 

y/  =  a  coefficient  fr'om  Wu  (1980) 

W  =  wind  speed 

&  =  wind  direction  counterclockwise  from  easterly 
h  =  depth 

=  x,y  velocity  components  at  the  water  surface 
4*  =  water  surface  elevation 
Ua,Va  =  x,y  velocity  at  the  bed 
a  =  bed  elevation 
s  =  salinity 

=  diffusion  coefficient  for  salt 
E  =  eddy  viscosity  components 
The  continuity  equation 


0 

dx  dy  dz 

is  solved  as  a  second  part  of  each  solution  step.  Equation  A6  is  converted  to  an 
appropriate  boundary  value  problem  through  differentiation  with  respect  to  z. 
After  rearrangement  it  takes  the  form 


d^'w_  d  f du  ^  dv^ 
dz^  dzydx  dy  ^ 


(A7) 


subject  to  boundary  conditions  specified  for  the  water  surface  and  the  bed. 


w/-  =  ur—+v^—+—  at  the  water  suface 
^  dx  ^  dy  dt 


(A8) 


and 
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da  da 

'Wa-Ua—+Va  — 
OX  oy 


at  the  bed 


(A9) 


Note  that  in  these  equations  the  values  of  m  and  v  will  be  known  at  all  locations 
from  the  previous  part  of  the  solution  step.  Values  of  w  in  this  solution  are  used 
in  the  next  iteration  for  u,v,h,  and  s. 

The  geometric  system  varies  with  time;  i.e.,  the  water  depth  h  varies  during 
the  simulation.  In  order  to  develop  an  Eulerian  form  for  the  solution,  it  is 
desirable  to  transform  this  system  to  one  that  can  be  described  with  a  constant 
geometric  structure.  Early  development  of  the  model  (King  1982)  used  a 
<T-transformation  in  which  the  bed  and  the  water  surface  are  transformed  to 
cons^ts.  In  a  later  analysis  of  this  method,  King  (1985)  pointed  out  that  at 
locations  where  a  sharp  break  in  bottom  profile  occurs,  the  transformation  is  not 
umque  and  momentum  in  the  component  directions  may  not  be  correctly 
preserved.  An  alternative  transformation  that  preserves  the  bottom  profile  as 
defined,  but  transforms  the  water  surface  to  a  constant  elevation  is  now  used  (z^ 
transformation). 

This  transformation  is  defined  by: 

^  (AlO) 

(All) 

z=a+(z-a)—^  (A12) 


where  b  is  the  fixed  vertical  location  to  which  the  water  surface  will  be  trans¬ 
formed.  Equations  A1-A6  and  A7-A9  then  incorporate  the  transformation  (AlO- 


Another  advantage  of  this  transformation  is  that  it  produces  =  constant 
lines  that  are  close  to  horizontal,  i.e.,  z  =  constant  lines.  This  results  in  less 
fictitious  density-driven  currents  near  bed  profile  breaks  (Stelling  and  van  Kester 
1993).  Since  stratification-related  phenomena  are  usually  nearly  horizontal,  it  is 
important  that  the  transformation  leave  constant  surfaces  that  are  nearly 
horizontal.  Considering  the  pressure  gradient  (due  to  the  density  gradient)  in  this 
transformation  produces 


dP_  dP  ,  dP  dz^ 
dx 


In  a  strongly  stratified  stagnant  system  this  pressure  gradient  should  be  zero. 
However,  note  that  Equation  A13  in  the  transformed  system  is  dependent  upon 
two  tenns  (each  of  which  could  be  large)  to  cancel  each  other.  This  could  cause 
artificial  currents  due  to  truncation  and  roundoff  error.  A  transformation  in 
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which  dz^  ldx«  a,  i.e.,  z''  «  z ,  will  reduce  this  problem.  Figure  A2  shows  an 
example  for  a  case  in  which  a  40-ft-deep  channel  passes  through  an  8-ft-deep 
bay.  Here  b  is  chosen  to  be  an  elevation  of  0  and  C  is  2  ft.  Near  the  break  in  the 
bed  profile  dz^  tdx  is  fairly  small,  or  z^  surfaces  are  nearly  horizontal. 
Contrast  this  with  the  a  transformation  in  Figure  A3.  The  a  =  constant  surfaces 
are  far  from  horizontal  along  the  channel  side  slopes.  The  truncation  and  roimd- 
off  errors  tend  to  drive  fictitious  currents  that  cause  the  denser  salt  water  to  leave 
the  channel.  The  z^  transformation  results  in 


ox 

whereas  the  a  transformation  is 


(A14) 


(A15) 


is  much  larger. 

The  Galerkin  finite  element  approximation  of  Equations  A1-A4  and  A7  uses 
a  quadratic  approximation  for  m,v,w,  and  s  and  linear  for  h  and  P.  The 
nonlinearity  is  addressed  by  Newton-Raphson  iteration  at  each  time-step. 
Generally  the  iteration  process  is  split  into  calculation  of  Equations  A1-A3,  then 
A7,  followed  by  A4.  This  sequence  is  repeated  until  sufficient  convergence  is 
reached. 
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Appendix  B 
Hurricane  Juan 


Introduction 

Based  on  the  results  of  the  initial  study,  additional  inquiry  was  requested  to 
evaluate  the  length  of  time  necessary  for  the  system’s  salinity  to  return  to  normal 
observed  levels  after  an  unusual  weather  event.  In  this  case,  the  surge  from 
Hinricane  Juan  is  used  to  constitute  this  event. 

To  accomplish  this  goal,  the  two  situations  previously  modeled,  open 
Mississippi  River  Gulf  Outlet  (MRGO)  and  closed  MRGO,  are  investigated 
again  under  extreme  conditions  by  applying  hurricane  surge  superimposed  on  the 
normal  astronomical  tide  to  the  easternmost  boundary  of  the  model,  in  the 
vicinity  of  Biloxi. 

Methodology 

The  tidal  record  for  the  NOAA  station  Cadet  Point  in  Biloxi  Bay,  MS,  was 
evaluated  to  pinpoint  the  passing  of  extreme  conditions  over  the  area,  such  as  a 
hurricane.  Hurricane  Juan  made  landfall  on  October  27, 1985,  and  was  used  for 
this  analysis.  The  observed  tidal  data  for  October  1985  were  then  compared  to 
the  predicted  tide  for  the  same  period.  The  predicted  astronomical  water  surface 
elevations  for  October  1985  were  generated.  Figure  B1  shows  both  the  predicted 
and  observed  water  surface  elevation  for  October  1 985. 

The  difference  between  the  two  data  sets  was  then  added  to  the  June 
astronomical  tide  used  in  the  prior  modeling  efforts.  The  resultmg  tidal  record 
was  then  applied  at  the  Biloxi  boundary  to  run  the  model  (See  Figure  B2). 

Five  locations  were  chosen  for  salinity  analysis;  the  middle  of  Lake  Borgne, 
Marcello  Castle,  Chef  Pass,  Little  Woods,  and  in  the  middle  of  Lake  Pontchar- 
train.  For  each  of  these  locations,  the  difference  between  normal  conditions  and 
the  hurricane  event  was  computed  under  both  situations:  open  MRGO  and  closed 
MRGO.  That  is,  the  salinity  for  normal  conditions  was  subtracted  from  the 
salinity  found  under  hurricane  conditions.  So  a  positive  salinity  means  that  the 
hurricane  raised  the  salinity  at  that  point.  This  was  done  for  MRGO  open  and 
also  with  MRGO  closed  and  are  shown  in  Figures  B3-B7. 
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Salinity  Difference,  ppt 


Figure  B5.  Difference  in  salinity  between  normal  and  hurricane  conditions  at  Chef  Pass 
with  MRGO  open  and  closed 


Figure  B6.  Difference  in  salinity  between  normal  and  hurricane  conditions  at 
Little  Woods  with  MRGO  open  and  closed 
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Mid  Lake  Pontchartrain 
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Figure  B7.  Difference  in  salinity  between  normal  and  hurricane  conditions  at  mid 
Lake  Pontchartrain  with  MRGO  open  and  closed 


Discussion 

The  figures  illustrate  that  the  area  near  the  MRGO  channel  are  more  affected 
by  the  hurricane.  Another  observation  is  that  there  is  generally  some  slight 
residual  of  the  hurricane  remaining  after  over  3,000  hr.  So  the  system  has  a  very 
long  residence  time.  The  difference  between  MRGO  open  and  closed  seems  to 
shrink  rapidly,  so  that  the  recovery  time  for  each  seems  to  be  about  the  same. 


Appendix  B  Hurricane  Juan 


- - - - - - -  T  Form  Approved 

REPORT  DOCUMENTATION  PAGE  _ | _ omb  No.  0704-0i88 

Public  reporting  bunten  for  this  collectpn  of  infoimation  is  |P  )  *1°", 

data  needed,  and  completing  and  reviewing  this  collection  of  infbrmatiOT.  noarations  and  Reoorts  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Ailingt^,  VA  222^^ 

valM  OMB  control  number.  PtJEASE  DO  NOT  RETURN  YOUR  F<^  TO  TOE  AB^ADDR^. - - -  3  DATbS  COVERED  ^ftom  -  TO;) 

1.  REPORT  DATE  (DD-MM-YYiV  2.  REPORT  TYPE  _ 

July  2001  _ I  Final  Rep^ - - - - - sa  CONTRACT  NUMBER 

4.  TITLE  AND  SUBTITLE  .  t,  _ 

Salinity  Changes  in  Pontchartrain  Basin  Estuaiy,  Loinsiana,  Resultmg  sb  GRANT  NUMBER 

from  Mississippi  River-GulfOutlet  Partial  Closure  Plans  _ _ _ _ 

5c.  PROGRAM  ELEMENT  NUMBbK 

_ _ _ _ _ _ _ _ — - 5d.  PROJECT  NUMBER 

6. AUTHOR(S)  ^ 

A.  R.  Carrillo,  R.C.  Berger,  M.S.SarrufitB.J.  Thibodeaux  se.  TASK  NUMBER - - - 

5f.  WORK  UNIT  NUMBER  “ 

7. Peero««INeORGAHUATK»IIM.E(^)«»««^  - TK^iiii^o^niSnioN^ORr- 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Engineer  Research  and  Development  Center 
Coastal  and  Hydraulics  Laboratory 
3909  Halls  Ferry  Road 
Vicksburg,  MS  39180-6199 


ERDC/CHLTR-01-14 


Vsp0e5orti«»MoMirote»e*GB.cYM«.bw««-u7UKh^ - io.sro««>»«oN,roK'SAcKu»m,b, 

U.S.  Army  Engineer  District,  New  Orleans 

P.O.  Box  60267  ii.  SPONSOR/MONITOR’S  repori 

New  Orleans,  LA  60160-0267  NUMBER(S) 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT  ~~ 

Approved  for  public  release;  distribution  is  unlimited. 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

The  Mississippi  River-Gulf  Outlet  (MRGO)  consists  of  a  ship  channel  36  ft  deep  and  500 
from  the  junctio^ofthe  Inner  Harbor  Navigation  Canal  and  the  Gulf  Intracoas^WWmNewCfr^^to^^^^^ 
waterVft  contour  in  the  Gulf  of  Mexico.  The  purpose  of  the  MRGO  is  to  provide  a  deep  draft  chaimel  to  &e  Port  of  New  OrleaiK  Im 
HaibOT  feciUties.  Since  the  MRGO’s  completion  in  January  1968,  saltwater  flux  from  the  MRGO  toou^  dir^ 

Borgne  and  the  Gulf  Intracoastal  Waterway  has  contributed  to  an  increase  in  &e  salinity  con^^on  ' 

This^port  presents  the  results  of  a  numerical  model  investigation  used  to  predict  av^e  sa^  c^ges  ^  ^ 

Lake  Pratctotrain  Basin  as  aresult  of  varying  levels  of  closure  of  the  Mississippi  River-Gulf  Outlet  below  Lake  Borgne. 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


15.  SUBJECT  TERMS 
Mississippi  River 

Gulf  Intracoastal  Waterway 

16.  SECURITY  CLASSIRCATION  OF: 


Gulf  outlet 
MRGO 


a.  REPORT 

UNCLASSIFIED 


b.  ABSTRACT 


c.  THIS  PAGE 

UNCLASSIFIED 


17,  UMITATION 

18.  NUMBER 

OF  ABSTRACT 

OF  PAGES 

44 

19b.  TELEPHONE  NUMBER  Cmcludearea 
code) 

Stancterd  Form  298  (Rev.  6-98) 
PnscrilMd  by  ANSI  Std.  239.18 


